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Abstract 

A periodic change of slow environmental parameters of a quantum system induces quantum 
holonomy. The phase holonomy is a well-known example. Another is a more exotic kind that 
exhibits eigenvalue and eigenspace holonomies. We introduce a theoretical formulation that de- 
scribes the phase and eigenspace holonomies on an equal footing. The key concept of the theory 
is a gauge connection for an ordered basis, which is conceptually distinct from Mead-Truhlar- 
Berry's connection and its Wilczek-Zee extension. A gauge invariant treatment of eigenspace 
holonomy based on Fujikawa's formalism is developed. Example of adiabatic quantum holon- 
omy, including the exotic kind with spectral degeneracy, are shown. 
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1. Introduction 

Consider a quantum system in a stationary state. Let us adiabatically change a parameter of 
the system along a closed path where the spectral degeneracy is assumed to be absent. We ask the 
destination of the state after a change of theparameter along the path. This question is frequently 
raised in discussions of the Berry phase [1]. An answer, which is widely shared since Berry's 
work, is that a discrepancy remains in the phase of the state vector, even after the dynamical 
phase is excluded. Indeed this is correct in a huge amount of examples |fl However, it is 
shown that this answer is not universal in a recent report of exotic anholonomies pi] in which the 
initial and the final states are orthogonal in spite of the absence of the spectral degeneracy. In 
other words, the eigenspace associated with the adiabatic cyclic evolution exhibits discrepancy, 
or anholonomy. Furthermore, the eigenspace discrepancy induces another discrepancy in the 
corresponding eigenenergy. 

For the phase discrepancy, an established interpretation in terms of differential geometry 
allows us to call it the phase holonomy |5]. This interpretation naturally invites its non-Abelian 
extension, which has been subsequently discovered by Wilczek and Zee in systems with spectral 
degeneracies Contrary to this, any successful association of the eigenspace discrepancy with 
the concept of holonomy has not been known. 
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The aim of this paper is to demonstrate that an interpretation of the eigenspace discrepancy in 
terms of holonomy is indeed possible. To achieve this, we introduce a framework that treats the 
phase and the eigenspace holonomies in a unified manner in Section|2] The key concept is a non- 
AbeHan gauge connection that is associated with a parameterized basis |01, and the identification 
of the place where the gauge connection resides in the time evolution. This is achieved through a 
fully gauge invariant extension of Fujikawa's formulation that has been introduced for the phase 
holonomy JMIlt]. Our approach is illustrated by the analysis of adiabatic quantum holonomies 
of three examples. First, Berry's Hamiltonian with spin-^ is revisited in Section|3] The role of 
parallel transport ifioi lsll. which accompanies the multiple-valuedness of a parameterized basis, 
in our formulation will be emphasized. The second example, shown in SectionlH exhibits exotic 
holonomies without spectral degeneracy. The last example, shown in Section |5] is the simplest 
examples of the exotic holonomies in the presence of degeneracy, i.e., the eigenspace holonomy 
a la Wilczek and Zee. Section |6] provides a summary and an outlook. A brief, partial report of 
the present result can be found in Ref. ifTHl . 



2. A gauge theory for a parameterized basis 

Two building blocks of our theory, a gauge connection that is associated with a parameterized 
basis 10], and Fujikawa formalism, originally conceived for the phase holonomy, are presented 
in order to introduce our approach to quantum holonomies. 

2.1. A gauge connection 

In the presence of the quantum holonomy, basis vectors are, in general, multiple-valued as 
functions of a parameter. In order to cope with such multiple-valuedness, we introduce a gauge 
connection for a parameterized basis. This has been introduced by Filipp and Sjoqvist |7] to 
examine Manini-Pistolesi off-diagonal geometric phase 11 1211 . As is explained below, this gauge 



connection is different from Mead-Truhlar-Berry's II13L [IJ] and Wilczek-Zee's gauge connec- 
tions i^, which describe solely the phase holonomy. 

For A^-dimensional Hilbert space 'K, let {\^n(s))]'^I^ be a complete orthogonal normalized 
system that is smoothly depends on a parameter s. The parametric dependence induces a gauge 
connection A{s), which isaN xN Hermite matrix and whose (n, m)-th element is 

OS 

By definition, A{s) is non-Abelian. For given A{s), the basis vector \^„{s)} obeys the following 
differential equation 

i-^lU^)) = y A„„(s)|^„(s)) (2) 
OS ^ 

n 

and we may solve this equation with an "initial condition" at i = i'. 

The dynamical variable of the equation of motion (|2]) is an ordered sequence of basis vectors, 
also called a frame, 

f{s) = [ms)\ 1^1 w), \^N-i{s))\. (3) 
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Its conjugation 



/(*)"'' = 



<^o(^)l 
(^i(^)l 



(4) 



is also useful. For example, the resolution of unity by {\^n{s))]^^Q is expressed as 

where l-^ is the identical operator for 'K, and the gauge connection A(s) is written as 

A(s)^i[f(s)^}-f(s). 



(5) 



(6) 



Now we have the equation of motion for f{s) 



i-f{s) = f(s)A(s). 
OS 



Its formal solution is 



/(/') = /(^')exp 



(-r 



A(s)ds 



(7) 



(8) 



where exp^ is the anti-ordered exponential for the contour integration by s Note that we 
need to specify the integration path to deal with the multiple-valuedness of f(s), in general. 

Our designation "gauge connection" for the whole A(s) is intended to clarify the difference 
from two famous gauge connections for the phase holonomy. When we choose {l^nis))]'^!^ as an 
adiabatic basis of a non-degenerate Hamiltonian H(s) with an adiabatic parameter s, the elements 
of the gauge connection A(s) have a well-known interpretation: A diagonal element A„„(i) is 
Mead-Truhlar-Berry's Abelian gauge connection for a single adiabatic state |^„(s)). The off- 
diagonal elements are nonadiabatic transition matrix elements, and, constitute "the field strength" 
corresponding to Mead-Truhlar-Berry's Abelian gauge connection [IJ. This also applies to a 
degenerate Hamiltonian with Wilczek-Zee's non- Abelian gauge connection, which describes the 
change in an eigenspace. On the other hand, A(s) that contains all elements {A,„„(.s)) is defined 
with respect to the change of the frame /(s), instead of each eigenspace |7]. As to be shown 
below, A(s) plays the central role in the quantum holonomy. 



2.2. An extended Fujikawa formalism 

Fujikawa has introduced a formulation to examine the quantum holonomy accompanying 
time evolution that involves a change of a parameter IMS- We will focus on the unitary time 
evolution for pure state in the following. As a building block of the time evolution, we exam- 
ine a parameterized quantum map, whose stroboscopic, unit time evolution from li/r') to is 
described by 

1^") = U(sM'}, (9) 

where U(s) is a unitary operator with a parameter s. This is because periodically driven systems, 
whose Floquet operator is U(s), are our primary examples, and our approach is immediately 
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applicable to a Hamiltonian time evolution, where U(s) correspond to an infinitesimal time evo- 
lution operator Let us vary s along a path C in the parameter space during L iterations of the 
quantum map (|9]l, where s' and s" is the initial and finial points, respectively. Accordingly we 
examine the whole time evolution operator 

f>({.^/)f=o ) = U{sL-i)U{sL-i) ■ ■ ■ U(so), (10) 

where si is the value of s at l-th step. Although the present formulation is applicable to investigate 
nonadiabatic settings, our primary interest here is an adiabatic behavior induced by the limiting 
procedure L — > oo. The "/(i)-representation" of the building block U(s) of the whole evolution 
isaN xN unitary matrix 

Z(s) = [f(sy)u(s)f(s). (11) 

In other words, we have U(s) - f{s)Z(s) {/(s) ' }. In order to deal with the change of s from 5/ to 
S/+1, we have 

U(si)=f(sM)[f(sM)^}xU(si) 

^f(s,M(sM,s,){f(s,y'), (12) 

where an effective time evolution matrix Zp(si+i, si) incorporates the unit dynamical evolution 
and the parametric change of s along a part of path C 



ZpisM,si) = sxp{i I Ais)dsjZ{si), (13) 
where exp^ is the ordered exponential. Hence the whole time evolution is expressed as 

u({s,}f:^) = f(s")Ba({si}lo) [fis'y} (14) 

where we have an "effective" time evolution operator 

fid({^;)io) = SL-i)Zp(sL-i,si-2) ■ ■ -Zpisuso). (15) 

Finally, we have "the /(i')-representation" of the whole time evolution operator 

U({s,}f-^) = f(s')W(C)Ba({s,}lo) {/(s')'} , (16) 

where 



H'(C) = exp|-/ J" A(i)t/ij. (17) 

Since the above definitions are exact, our formulation is invariant against a basis transforma- 
tion with N X N unitary matrix G{s) 

f(s) ^ /(i)G(i) (18) 
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once we incorporate the following transformations 



A(s) ^ G(s)^A(s)G(s) + iGisY^^^, (19a) 

OS 

W(C) ^ G(s'yW(C)G(s"), (19b) 

BMsi}io) ^ G(s"yBMs,}io)G(s'). (19c) 

This is Fujikawa's hidden local gauge invariance |9] in a generalized form. The strategy of 
Fujikawa formalism is to extract a geometric information from the whole time evolution operator 
via the expression ( fT6t with an appropriate restriction of G(s), as is shown below. 

Let us examine the case that the change of s is slow enough so that we may employ the 
adiabatic approximation |fl4ll . Accordingly it is suitable to choose f(s) as an adiabatic basis, 
i.e., each basis vector \^„{s)} is an eigenvector of U(s), to make Z(s) a diagonal matrix, whose 
non-zero elements are the eigenvalues of U(s). Let Zn(s) be the eigenvalue of U (s) corresponding 
to an eigenvector \^„(s)), i.e., 

UisMn(s)) = Z„(sMn(s)}, (20) 

where we assume that there is no degeneracy in eigenvalue. Note that Zn(s) is unimodular due 
to the unitarity of U(s). Now the gauge transformation G(s) is restricted to U{1)'^^ times a 
permutation matrix, which correspond to the freedoms to choose the phases of basis vectors, and 
assign the quantum numbers, respectively. The permutation matrix is required to deal with the 
eigenspace holonomy, as is shown below. In terms of Zpisi+i, si), the adiabatic approximation is 
the diagonal approximation dgj 

ZF(i,+i,i,) -Z°(^,+i,i,) = expJ/ J^"'A°(5)flfiJz(i,), (21) 

where A°(.s) is the diagonal part of the gauge connection A(.s), i.e., A)J5„(s) - 6mnAmm(s)- Namely, 
A^„,(s) is Mead-Truhlar-Berry's gauge connection for the m-th state \^m(s)} 1 1 3l l If . The corre- 
sponding adiabatic approximation of £d({'S/)fLo) is 

B-A{si]io) = Z^iSL, SL-l)Z^iSL-USL-2) ■ ■ .Z^(suSo), (22) 

which are decomposed into two parts: 

fiad({*/)io) = B(C)D({s,}l,,), (23a) 

where 

B(C) = exp^/J" A°(i)c/ij (23b) 

is the geometric part, and, 

L-l 

D({s,}l,) = Y]z(si) (23c) 
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contains dynamical phases. We retain the path-ordered exponential in B(C) to make it applicable 
to the cases with the presence of spectral degeneracies, as shown below. To sum up, the adiabatic 
approximation of the whole time evolution is 



(24) 



where 



M(C) = W(C)B(C) 

- exp j' A(s)d. 




Jc 



(25) 



is the geometric part determined by the path C, and two gauge connections A(s) and A (s). 

For a cyclic path C, we call M(C) (Eq. (l25l l) a holonomy matrix, which describes the adiabatic 
change of state vector, starting from an eigenstate at s = s', along the closed path C. 

An explanation why W(C) is required to describe the eigenspace holonomy is the following. 
Let us assume that /(s) is single-valued. This implies that WiC) is the N x N identical matrix. 
Consequently, M(C) = B{C) is always diagonal and thus cannot describe the eigenspace holon- 
omy, though the single-valuedness assumption does not prevent the conventional approach from 
describing the phase holonomy. On the other hand, in the presence of the eigenspace holonomy, 
a factor of M(C) need to be a permutation matrix. Since B{C), which is always a diagonal matrix 
according to its definition, cannot be such a factor, the permutation matrix is need to be supplied 
by W(C). This is consistent with the fact that the presence of the eigenspace anholonomy implies 
the multiple-valuedness of f(s). 

Furthermore, when we employ the parallel transport condition A^{s) — ifioi [sl 01, the 
holonomy matrix takes extremely simple form: 



which is determined only by W(C). In other words, all the adiabatic quantum holonomies can 
be summarized as a holonomy in the orderd basis f(s) (Eq. (O). For the phase holonomy, this 
observation is akeady reported by Fujikawa fllst]. In this sense, the parallel transport condition 
offers a privileged gauge. 

We explain the consequence of the gauge transformation (fTsT l under the adiabatic approxi- 
mation, where G{s) is restricted to be a product of a permutation matrix and a diagonal unitary 
matrix. The invariance of the adiabatic time evolution operator, which appears at the right hand 
side of Eq. (|24] |. is assured due to the following 




(26) 



A°(s) ^ G(syA°(s)G(s) + iG(s) 



dG(s) 
ds 



(27a) 



WiC) ^ Gis')^W{C)G(s"), 
B(C) ^ G(s"y B(C)G(s'). 



(27b) 
(27c) 



Hence we obtain the manifest covariance of M{C) 



M(C) ^ G(s')^MiC)G(s'). 



(28) 
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Eq. i25[ provides a correct expression of the holonomy matrix M(C) not only for maps ^ 
but also for flows, i.e., Hamiltonian systems and periodically driven systems. For a system 
whose time evolution is generated by a nearly static Hamiltonian H(s), we will derive Eq. (IZST i 
in Appendix lAl with a suitable discretization of time. For a periodically driven system described 
by a Hamiltonian H{t, s), where H{t, s) = H(t + T, s) is assumed, a Floquet operator 




is the unitary operator U(s) to provide a stroboscopic description of the system. Hence this 
system is reduced to a quantum map. 

An extension of our formulation to the case that the presence of spectrum degeneracy whose 
degree is independent with s along a closed path C is shown. The resultant expression for the 
holonomy matrix dZST l remains the same. This is achieved by a suitable extension of gauge 
connections A{s) and A^{s). For the eigenspace corresponding to the eigenvalue z„(i) of U{s), 
we have a normalized orthogonal vectors \^ny(s)) with an index v for the eigenspace, where 

UisM„y(s)) = z„(sM„y(s)} (30) 

and {^n"v"(^Mn'v'is)} - 6„"„'6y"y'. The gauge connection for the parameterized basis is 

d 

An"y",„'y'(s) = i{^„"y"(s)\ — \^,fy{s)). (31) 

OS 

It is straightforward to see that A^{s) that appears in Eq. dZSl l is non-Abelian: 

■^«"v",«'v'(*) — ^n"n'An"v" ,n'v'i.S), (32) 

where Aj^^„ „y,{s) is Wilcek-Zee's gauge connection for the n-th eigenspace |0]. 



3. An example of the phase holonomy 

Our formulation is applied to Berry's simplest example of the adiabatic phase holonomy fill- 
Let us suppose that a spin-^ is under static magnetic field B. With a suitable choice of units, the 
spin is described by a Hamiltonian 

H{B) &, (33) 

where & = 2 j=.v,y,z j is the Pauli operator for the spin, and, ej (j = x, y, z) is the unit vector for 
7-axis. The spectrum of H{B) is {+B}, where B = \\B\\. To investigate the adiabatic holonomy, we 
have to exclude the degeneracy point B = 0. The unit vector n = B/B is parameterized with the 
spherical coordinate, i.e., n = cos (fsm6 + e,, sin (psinO + e, cos 0. Let \^±{B)) be a normalized 
eigenvector of H{B), corresponding to the eigenvalue +B: 

\UB)} = e--'^'^ cos ^1 T ) + e''^'^ sin ^| i ) (34a) 
\UB)} = -e-''f'^ sin ^1 T ) + e''^'^ cos ^| i ). (34b) 
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Note that \^±(B)) are multiple-valued as functions of (6,(p). To assure them single-valued, the 
range of (0, ip) needs to be restricted within an open set Um = {{S, (p)\0 < 6 <n and < (/? < 2n], 
for example. The corresponding frame f{B) = is 



/(B) = [lT>, li>l 



e-'fl^ cos f 
e'^l^ sin f 



-e-'fl^ sin f 
e'fl^ cos f " 



(35) 



Gauge connections Aj^(B) = i\f{B)^\df{B)ldx{x - 0, i^, B) for parametric changes of B are 



Ae(B) = -a,, 

Aip(B) - - {cr~ cos - cTf sin 0) , 



Ab{B) = 0, 

where we employ 2x2 complex matrices 








1 




-i 




1 


CTx = 


1 





, cr,. = 


i 


cr, = 


-1 



Accordingly, the Mead-Truhlar-Berry gauge connections are 

A^(B) = 0, A°(B) = i(r,cos0, AD(B) = 0. 



(36a) 

(36b) 
(36c) 



(37) 



(38) 



As is well known, the strength of the magnetic field B plays no particular role in the computation 
of the phase holonomy once B is kept nonzero. 

Evaluations of the holonomy matrix M{C) for typical closed loops in the parameter space are 
shown. First, we examine a loop C in which f{B) is single-valued, e.g. C c Um- Consequently 
WiC) is the 2x2 identical matrix. This is a conventional wisdom to obtain a formula of the 
phase holonomy, where an appropriate gauge for f{B) (or, equivalently, \^±{B))) against the 
loop is chosen to avoid the multiple-valuedness of f{B) [1]. Hence all the holonomies reside in 
B{C) - M(C). The evaluation of the contour integral in B(C) is straightforward to obtain the 
classic result in the matrix M(C): 



M(C) = exp 



(-iQ(C)tr,), 



(39) 



where f2(C) is the solid angle for C iH. 

Second, we examine a meridian great circle Cg in which moves to In with ip kept fixed. 
The parametric change along such a circle can induce a change of the sign of an eigenvector Ml 511 . 
In our example, f(B) cannot be single-valued on Cg: 



f(B)\ 



f(B)\ 



0]2k ■ 



(40) 



Although the conventional strategy mentioned above is to avoid such multi-valuedness, we insist 
the present choice of the gauge ( l35T l to show an alternative way to reproduce the conventional 
result. Thanks to the present choice of the gauge, the Mead-Truhlar-Berry gauge connection 
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satisfies the parallel transport condition - and B(Ce) - 1. This enable us to employ 
Eq. (l26T l to obtain the holonomy matrix: 

M(Cfl) = exp{-i j) ^-(Tyde^ = e-""^' = -1, (41) 

which is consistent with Eq. ( [39] l. This is an example that our formulation properly deals the 
multiple-valuedness of f{B). On the other hand, if we choose an appropriate gauge to make 
f{B) single-valued on the circuit Q, W{Ce) is trivial to put all the nontrivial holonomy in B(Cg), 
as is stated above. 

Finally, let us consider a circle of latitude C^, where is kept fixed and (p is increased from 
to Itt. It is straightforward to obtain B(C^) - exp( otct, cos 0) and W(C^) - -1. The latter 
indicates again the sign change in the parametric dependence along C^. These elements are 
combined to reproduce a well-known result 

M(C^) = exp {-in{\ - cos 0)a-^} . (42) 

We remark that, in all the examples above, the holonomy matrices M are diagonal, so that the 
eigenspace holonomy is absent. Accordingly the eigenvalue anholonomy is also absent. 



4. An example of the exotic holonomies in quantum map spin-| 

In Berry's Hamiltonian (Eq. (l33Tl). the strength of the magnetic field B plays no role in quan- 
tum holonomies. One reason is that the corresponding gauge connection A ^(B) vanishes. An- 
other reason is that it is impossible to make any loop in the parameter space by an increment 
of B, with n being kept fixed. To make a loop for a strength parameter, we may examine the 
following quantum map for a spin-j 

exp {-iA (a + bn ■ &)} , (43) 

where « is a normalized real vector, a and b are real constants to be specified later, and A is 
the strength. The periodicity of the quantum map with respect to the increment of A implies 
that there is a loop in the parameter space of A. In particular, if we choose a - q/2 and b - 
(2 - q/2) with an integer q, Eq. (l43T l is periodic as a function of A, with a primitive period 
2jt. Accordingly, the parameter space of A is identified with S ' and it might be suitable to 
investigate quantum holonomies for a periodic variation of A. However, such a loop does not 
allow us to study adiabatic holonomies, since there remains a spectral degeneracy along the loop 
at A - (mod 2n). 

A simple way to lift the degeneracy at /i = is to concatenate two quantum maps: 

exp ^-ifi + ■ o-jj exp |-/^ J| + ■ o-J| , (44) 

where q and p are integers, m and n are normalized vectors in M^, and, yU and A are strengths. Due 
to the periodicity in yU and A, the parameter space of (/i. A) is a two-dimensional torus 5 ' x 5 ' . 
Both m and n specify points on a sphere S^. In the following, we fix m = e, and parameterize n 
by spherical variables and ^ as n = c i cos tp sin 9 + e,. sin (fsm9 + cos 0. If we change m with 
keeping m ■ n fixed, it induces only the Berry phase. 
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To facilitate the following analysis, we examine the symmetric version of the quantum map (l44l l 
U . exp {-,1 exp l^-iA [^ + ^^n-cr 

xexp|-/-(- + ^^/M-o-jl. (45) 



2\2 2 

For brevity, we omit the parameters in the following. A possible implementation of the quan- 
tum map i45[ is available by a periodically driven system that is described by the following 
Hamiltonian 

^(0 = J" (f + ■ + ^ (f + ^« ■ ^] Z - ^■)' (46) 

where the Floquet operator for a unit time interval -1/2 < f < 1/2 is f). The magnitudes of the 
magnetic fields of the unperturbed system and the perturbation are 



B^, = ^i2- q)n, 
Ba = \{2 - p)A, 



(47) 



respectively. It is straightforward to show 



U = e-'^^^^P^'^ Jcos ^ - /o- ■ /j , 



(48) 



where 



A s 2 cos ' (cos Bfj cos - cos 6 sin sin B^j , (49a) 
/ = (sin Bfj cos Ba + cos cos Z?^ sin Z?,i) /w 

+ (sin Ba) [n-{n- m)m] . (49b) 



It is also easy to see - sin^(A/2). Hence / = // sin(A/2) is a unit vector Now we have 



and its eigenvalues are 



^/.exp^,-(^.f^./|^ (50) 



Z..exp^/(^4]r (51) 



Corresponding quasienergies are E± = {/iq + Ap + A)/2, which is defined up to modulus In. 

In order to study the adiabatic holonomies, we need to identify the spectral degeneracies, 
whose condition is e"^ = 1, in the parameter space. It is useful to see A as a function of Ba 

A = 2 cos"' (A cos (Ba + a)) (52) 
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where A = ^1 - sin^ 6 sin^ and a is an "initial phase" that is independent with A. We choose 

the branch of cos"' A in [0,7r]. If A < 1 , A oscillates within the range [2cos ' A,27r-2cos~' A] c 
(0, In) as a function of B^, and encounters no spectral degeneracy. On the other hand, the con- 
dition A = 1 (i.e., sin sin = 0) implies the presence of the spectral degeneracy. From the 
similar argument for B/j, we will encounter spectral degeneracies if sin sin sin = 0. 
Let us examine the case sin = 0. Since this implies cos = +1, we have 

A = 2 cos"' (cos [b^ ± B^)) . (53) 

Accordingly the degeneracy points draw lines in (B^, B^)-plane as 

iB^,Bd\— (54) 



corresponding to the condition cos 6 = ±1. On the other hand, if we assume sin B^ sin B^ = 0, 
we have 

A = 2 cos"' (cos Bjj cos B^) . (55) 

Hence another condition for the spectral degeneracy is | cosB^ cosB^I = 1, i.e., the degeneracy 
points are at lattice points: 

(B^,BM^,^el\, (56) 
^ 'inn ) 

for all 6. 

Summarizing above, we show the location of the spectral degeneracies in terms of (jj., v), 
whose space is the two-dimensional torus. The degeneracy lines are specified by (ji. A, 9) as 

Ou, A, 61) G Z, - e Z ^ . (57) 

I In ^ ) 

In addition to this, we have isoleted degeneracy points as 

\^~.^^~\Wl-A- (58) 



2-q 2-p 

Except these degeneracy points, it is legitimate to introduce a zenith angle of /, s.t., 

sin B„ cos B^ -I- cos cos B„ sin B^ 

cos© = 



sinO 



sin(A/2) 
sin sin B^ 



(59) 



sin(A/2) 

It is straightforward to obtain the eigenvectors of U, corresponding to the eigenvalues z±: 

1^^) = e"'^/2 cos 1 1 T ) + sin 1 1 i >, (60a) 

|^_) = -e-'-''^ sin || T > + e'^'^ cos |u ). (60b) 
11 



Let / = [|^+), be a frame. The gauge connections Q for / are 

1 50 



Ag 

A 



- (0-7 cos - cTj. sin 0) , 



1 50 

2 ■ dA 
1 50 

^'^ ^ 2"> ' 



(61a) 
(61b) 
(61c) 
(61d) 



and the corresponding Mead-Truhlar-Berry gauge connections are Ag - 0, - jcr^ cos 0, 
A° = 0, and - 0. With these gauge connections, we will examine the holonomy matrices of 
typical loops in the parameter space. 

First, we examine the meridian great circle Cg in which 6 moves to In with other parameters 
are kept fixed. It is straightforward to see B{Cg) = 1, due to the parallel transport condition 
Ag - 0. Hence all the holonomies reside in W(Cg) = M(Cg): 



1 



= exp(-/-cr„0|-o|, 



(62) 



where @\g1o^ the change of along Cg, is determined by the image of Cg in the sphere (0, ip). 
If the image is a closed circle, we have ©IIZq - ±2;r, where + correspond to the direction of the 
path. Both cases provides a Longuet-Higgins type phase change M(Cg) - e™"''' = - 1 . On the 
other hand, if the image is closed self-retracing curve along an ark, we have 0|g^o ~ 0' which 
implies M{Cg) - 1 . The following index 



(63) 



Ba + B^ 




Ba - B^ 


n 




n 



where [x\ is the maximum integer not greater than x, determines which is the case, as shown in 
AppendixiBl 



1 2n{-\)''^^ for r is even 



for r is odd 



(64) 



Hence we obtain 



M(Ce) = (-1)' 



(65) 



Next, we examine a circle of latitude C^, where ^p is increased from to 2n and the other 
parameters are kept fixed. We have B(C^) - ex-]){incr^ cos 0) and W(C^) = -1. Accordingly, we 
obtain 



M{C^) - exp (-iniX - cos 0)crj . 



(66) 
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So far, the holonomy matrices M(C) are diagonal, so that neither Cg nor incorporates the 
exotic holonomies. The following is the first example of the exotic holonomies in this paper. 

Let us examine a closed loop Ca, in which A is increased from to 2n, being kept fixed other 
parameters. To avoid degeneracies along Ca, we choose fi{2 - q)/(27T) i Z and 0Iji i Z. When 
we increase A from A - A' io A - A' + 2n, we have 

A I /'^l'i=''' for even /7 

I ZTT - A|^^^/ for odd p 

i.e., an anholonomy in A occurs. Accordingly we have an eigenvalue holonomy 

I {z±\a=a' for even/? 

\z^\a=a' for odd/? 

which implies the presence of the eigenspace holonomy. We proceed to evaluate the holonomy 
matrix M{Ca)- Because of the parallel transport condition = 0, we have M{Ca) - W(Ci). On 
the other hand, we have 



IV(Ci) =exp(-/(0 ^(^y^dA 



•X 1 ' 



= exp|-ii£r,ei;j:J"). (69) 



Hence we need to examine the zenith angle 0. From Eg. (159b and sin(A/2)| j_j,+2n^ = sin(A/2)|^=^/, 
we have e'®\^^^^ - ^'*^L=o' ®li=o - - p) mod 2n. Since this is not suffice 

to determine the precise value of W{Ca), we need keep track of against the increment of A 
to obtain 0|^^o. Now the parameter space perpendicular to Ci is (0, i^, //) G 5^x5', which is 
divided into subspaces by the spectral degeneracies p - 2nkl{2 - q) (k is integer). Since, in each 
subspace, ©{^Zo is constant, it is suffice to evaluate it at a representative point. Let us choose a 
point 9 = 7r/2 and p - n(2k + l)/(2 - q), where the spectral gap takes a constant value A - n, 
from Eq. ( 149 at . Accordingly we have e'® = (-1)* exp {/(-1)*^(2 - p)/2] , from Eq. Hence 
we obtain 

0|5:o = (-1)^2-/7), (70) 

which also holds for 2nkl{2 - q) < p < 2n{k + l)/(2 - q), i.e., k - {p{2 - q)/(2n)]. Hence we 
have 

M(Ca) = exp J-/^(-l)V2 - p)cr, 

-cos^, -(-1)* sinful 
(-l)*sinf, -cosf J- ^'^^ 

In particular, if p is odd, the off-diagonal elements of M{Ca) remains, so that the eigenspace 
holonomy exhibits. This is consistent with the emergence of the eigenvalue holonomy for odd p. 

Now the similar analysis of quantum holonomies for the circuit C^, where p is increased from 
to 27T, is trivial. Hence we show only the holonomy matrix 



M(C„) 



-cosf^ -(-l)*^sinY 
(-l)*sin^, -cos^ 



(72) 



where k - [A(2 - p)l{2n)\. We conclude that odd q along implies the exotic holonomies. 
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5. Example 3: the exotic holonomies a la Wilczek-Zee 



A simple example of the eigenspace holonomy accompanying spectral degeneracy is shown. 
Extending Mead's study |3 12\ on non-Abelian adiabatic phase holonomy [6], we introduce a 
quantum map with Kramer's degeneracy. 



5. 7. Quantum map for spin- j with Kramers ' degeneracy 

To introduce our model, we review the time-reversal invariance structure in an atom with 
odd-number of electrons |.16, 17]. For a comprehensive explanation, we refer Avron et al. iHtIi . 
Let / be the total angular momentum of our system and \ J, M) the standard basis vector for /, i.e., 
f\J, M) = J(J+l)\J, M), J-e~\J, M) = M\J, M>, and J -(e.+ie,)]], M) = ylJ(J +\)-M(M± \)\J, M+ 
1). The standard time-reversal operator for / is an anti-unitary operator K = exp{-i7rJy)K(), 
where Kq is the complex conjugate operation in the |7, M)-representation. We examine the 
fermion case - -1, which implies that 7 is a half-integer If an Hermite operator commutes 
with K, its spectrum exhibits Kramer's degeneracy. The same is true for unitary operators. 

We focus on the case J - ^, and introduce basis vectors as follows: 



3 3 



\e2} = 



1, 



\Kei} = K(\ei})^\^,-^}, 
\Ke2} = K(\e2})^\l,h- 



(73) 



'2 2' ' ' " " '2 2' 
Our physical observables are spanned by the following time-reversal invariant operators 



TO = kiXeil + \Kei}{Kei\ - \e2){e2\ - \Ke2){Ke2l 

T, = \ex){Ke2\ + \Ke,){-\){e2\+h.c., 

T2 = \ei)i-i){Ke2\ + \Kei)(-i){e2\ + h.c, 

T3 = \ei}{e2\ + \Kei}{Ke2\ + h.c, 

t4 = |ei)(-0<e2l + \Kei}i{Ke2\ + h.c, 



(74) 



which are traceless and form a Clifford algebra 



Several properties of t„ are shown in AppendixICl 

We introduce an extension of the quantum map for spin-^ (Eq. ( |45]) ) 



t/ = exp ^ -/^ 1 1 -1- ^Y^To I exp <; -iA 



.U Iq 2 - <7 , 



P_ 2- j? 
2 2 



a=0 J) 



(75) 



(76) 



where (na,)*^^ is a unit vector in R^, i.e., 2„ = 1, and (q, p) e I?. The quantum map ( l76l l can 
be implemented by a periodically pulsed driven system in a similar way shown in the previous 
section for the quantum map (05]). Since the unitary operator ( l76b is 27r-periodic both in ji and 
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A, the parameter space of (fj.. A) forms a two-dimensional torus 5 ' x 5 The unit vector (n„) is 
parameterized by spherical variables: 



no - cos ( 
In AppendixiDl we show 







cosx 




«3 




COS (fi 




sin 77 sin 0, 




"2 




sinx 


«4 




simp 



cos // sin 6. 



U = e 



where A is defined in Eq. ( |49a| ). 



Iq - sin Bp cos B \ + cos cos Z?^ sin B,\, 



(77) 



(78) 



(79a) 



and 



la = ria sin 



(79b) 



for a ^ 0, where the definitions of B^ and are shown in Eq. ( |47] ). Since 2^^o 'a - sin(A/2)^, 
we normalize /„: 



1 



sin(A/2) 



la. 



Accordingly we have 



U - exp I -i 



+ 2 Zj 



(80) 



(81) 



The eigenvalues of U are 



z± = exp -I 



f^q + Ap ^ A 

2 "2 



(82) 



Namely, the spectrum is completely same with the example shown in Section HI To examine 
eigenvectors, we parametrize l^, with the zenith angle 0, which is aheady introduced for the 
quantum map spin-i in Eq. ( |59] l. We have 



/o - cos ©, 



From Appendix|El the eigenvectors of U are 



/l 






sin J] sin 0, 


h 




cos If 


h 




h 










sin(f 



cos //sin©. 





© 





|t/i)cos — + 


|d/2)sm-, 




© 


© 




\Kdi)cQS — 


+ \Kd2) sin — , 




1 © 


\ 


l'^i>|-sm- 


l + |t/2>COS-, 






0\ 




\Kdi)^ cos 


- + \Kd2) cos 



© 
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(83) 



(84) 



where 



(85) 



\di) = cos I) + \Ke,) sin |) , 

|c/2> = |e2> (e'f'^-^'/^ cos |) + \Ke2) (e-'(^--»^>/2 ^j^^ |J ^ 

ITTt/i) = |ei> gijj |j ^ l/Tei) cos ^) , 

\Kd2) = |e2) (-e+'f-^-^'^^ gjjj + i^g^^ {^g-i(f-x)a cos . 

Note that we put each basis vector before its complex coefficient above to prevent a confusion 
due to the presence of anti-Hermite operation K. To conclude this subsection, we introduce a 
frame composed by the eigenvectors tj ifisll : 

/=[|^+>, \K^^\ |^_>, \K^-)\. (86) 



5.2. Analysis of adiabatic holonomies 

We examine the adiabatic holonomies of the quantum map. Note that A and depend only 
on ji. A, 9 and is the same ones for the quantum map spin- 1/2 (see, Eqs. ( l45T l. and (|59l)). Hence the 
degeneracy points in the parameter space and the holonomy in the eigenvalues are the completely 
the same. We will focus on the eigenspace holonomy in the following. 

It is straightforward to obtain the gauge connection from the eigenvectors. Since © depends 
on ji. A, 0, the corresponding gauge connections are defined through the derivative of / by 0: 




+il2 



-ih 




5© 

dx ' 



(87) 



for X - 9,11, A, where I2 is the 2x2 unit matrix. The corresponding Wilczek-Zee gauge connec- 
tions vanish to satisfy the parallel transport condition, i.e.. A° = A° = A° = 0. For other gauge 
connections, ju. A, and 0-dependences are introduced through 0: 



A^p — — 



-cTy COS cr,. sin 

(Ty sin © CTy cos 

cr, cos -cr- sin 

- (T, sin -cr, cos 



1 

cos 77 H 

' 2 







1 

cos n -I — 
' 2 



CTx COS 

- cTv sin 









-cr^ sin 
-CTv cos 



sm 77, 



sm T], 



(88a) 
(88b) 
(88c) 



and 



A" - - 
^" " 2 



4° 



A" - - 



-cr, 








O-y 









COS 0, 



COS COS?; + 



1 



cos 77-1 

' 2 



cr., 
-CT;, 





a-x 



sm 77, 



cos sin 77. 



(89a) 
(89b) 
(89c) 
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It is straightforward to evaluate the holonomy matrix M(C„) JZSl l for a closed loop C„ where 
a(=fj., A, 0, T], (p,x) is increased from to In, once we take care the "anholonomy" in along C„, 
which is also clarified in Section|4] For jj. A, 6, we have 








-ih 


H 


+il2 






h COS 

+I2 sin 



2 



-/2 sin ^ 

/2C0S^ 



where 



01 



2k 
^=0 



{2n{- 1 )''^^ for r is even 
for r is odd 

(-l)[^(2-?)/(2;r)]^(2-;,), 



0|2^Q = (-1)[''(2-P)/(2t)]^(2-^), 



(90) 



(91a) 

(91b) 
(91c) 



and r is defined in Eq. ( l63T l. Hence M{Cg) exhibits only Herzberg and Longuet-Higgins' sign 
change BlSll 



M(Q) = (-l)'+74. 



(92) 



Also, the same kind of sign change appears along (Ci) with even p (q), i.e., M(C^) - 
(_l)i+p/2/4 (M(G) = (-l)'+«/2/4). On the other hand, a mixture of the eigenspace holonomy 
and the Herzberg and Longuet-Higgins' sign change occurs along with odd p 



M(C,) = (-1) 



W2-q)/(2K)]+(p-l)/2 



and, along Ca with odd q 



M{C,) = (-1) 



_-,x[,)(2-p)/(2;r)]+(?-l)/2 



-h 
h 



-h 
h 



(93) 



(94) 



which do not incorporate mixing within the degenerate eigenspaces. Other holonomy matrices 
M{Ca) (a = T], tp,x) describes genuine Wilczek-Zee's phase holonomies: 



exp (/(r,Q,/2) 

exp(-/cryQ;,/2) 

exp [-/(cr, cos rji + cr^ sin ?7i)Qi /2] 

exp [iXcTj. cos rji - (Tx sin ?7i)Qi /2] 

exp [-/(cr. cos 772 + (Tj, sin 772)02/2] 

exp [-/((T7 cos 7/2 - cTv sin 7/2)^22/2] 



M(C,) 
M(C^) 

M(C;,) 

where Q.,j 

P\ = a/1 - (sin cos 77)2 and yS2 = - (sin sin 77)2. 
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(95a) 
(95b) 
(95c) 



= 2;7r(l - COS0), = 2;7r(l - fij), rjj = -/ In |(cos cos 7; + / sin7;)/j6j} for y = 1,2, 



6. Summary and outlook 



We have introduced a framework that is capable of describing the eigenspace holonomy and 
the phase holonomy in a unified manner. Several examples have been shown. 

In hindsight, it might seem rather odd that, in two decades since the first discovery of Berry 
phase, the full, gauge invariant formulation for quantum holonomy has not been conceived prior 
to this work. It should probably be attributed to the lack of the incentive to improve on the 
original expression of Berry; Nothing other than the phase anholonomy has been anticipated for 
the adiabatic cyclic variation of parameters for regular Hamiltonian system. As a result, the gauge 
invariant formulation must have seemed a redundant luxury. However, if we once recognize the 
possibility of eigenstates' exchange without level crossing for cyclic parameter variation, both 
for singular systems and for time-periodic systems, it becomes imperative to treat the choice 
of basis frame explicitly within the formalism, which has naturally lead us to arrive at the fuU, 
gauge invariant formulae. 

In a simple minded view, it is the periodicity of quasienergy, which is a result of the time- 
periodicity of the system, that enables the eigenvalue holonomy in a natural manner For a 
Hamiltonian system, with the energy defined on entire real number, an energy eigenstate after 
cyclic variation of parameter cannot reach another eigenstate of different energy in a usual way, 
since the crossing of levels is prohibited for adiabatic variation. The only possible exception 
appears to be singular Hamiltonian systems, for which the highest and the ground eigenenergy 
diverge [4]. 

In this work, we have focused on the most elementary setting of quantum holonomy, i.e., the 
adiabatic excursion of pure quantum eigenstates along a closed path in the parameter space. A 
vast amount of studies on the phase holonomy naturally suggests possible directions of extension 
of the present result. We mention only few of them. A straightforward extension is to examine 
noncyclic path [ 19lfl2| . Also, loosening of the assumption of adiabaticity and resulting extension 
into e.g. Aharonov-Anandan's nonadiabatic settings 1I20I1 . are expected to be straightforward 
thanks to the generality of Fujikawa's formulation [9], which is the basis of our theory. It seems 
timely as well as interesting to examine the eigenspace holonomy in dissipative systems, for 
which we will need to horn appropriate techniques to treat of the eigenspace holonomy in mixed 
states [21]. 

Finally, we mention a question that is raised from the main result Eq. dZSb . which supplies a 
complete prescription to quantify the adiabatic quantum holonomy. How this helps to understand 
the exotic holonomies intuitively? Is it possible to find any underlying object or concept that 
governs them, for example, in the manner of diabolic point for the case of Berry phase? Herzberg 
and Longuet-Higgins has shown that the phase holonomy along a closed loop C implies the 
presence of spectral degeneracy in a surface S enclosed by C in the parameter space |15] (see 



also, Refs. I10ll22ll ). Is there any counterpart of the argument of Herzberg and Longuet-Higgins 



for the exotic holonomies? This does not seem likely, for the case of exotic holonomy, at the first 
glance, since there is no room to make S from C in all the examples shown in this paper We 
now believe that an affirmative answer is to be found in the complex parameter plane, on which 
we shall focus our attention in a forthcoming publication. 
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A. The gauge theory for Hamiltonian tune evolution 

We will derive the holonomy matrix (Eq. dZSt ) for a system described by Hamiltonian H(s) 
that depends on a time-dependent parameter s. Suppose that s is moved from s' to s" along a 
path C, during < t < T. The corresponding time evolution operator is 

f>({«rW]) = expj-i j| H(s,)dty (96) 

We divide the time interval [0, T] into L parts. Let f/ = {l/L)T. A short time evolution operator 
Ui is accordingly introduced: 

Ui = exp{-^jJ H(s,)dty (97) 

Wlien L is large enough, we have 

f), = exp(--^^(s,)e) + (9(e2), (98) 
where si = 5(5,,^, + s,,) and e = T/L. Now we introduce a i-dependent unitary operator 

f/, = exp(-^^(s)e). (99) 

Accordingly we have 

U({s,} ie[0,T] 

)-Y\ ^^.+0{e). (100) 



/=o 



Since we choose L so as to satisfy e <«; 1, our formulation explained in Sectionals straightfor- 
wardly applicable. 

First, for a given /(.?), the /(i)-representation of U(s) is 

Z(s) = exp(--^ {/(*)■'■} H(i)/(.)e) +0(62). (jQj) 
On the other hand, we have 

exp [i A(s)ds'^ = exp (iA(s,)se) + O (e^) , (102) 
where we assumed s = (9(1). Hence we have 

exp^/ A{s)d.^Z(s,) = exp(-^F(s,, s)e) +0[€^) , (103) 

where we introduce Fujikawa's Hamiltonian matrix ^ 

F(s, i) = {/(.) l } ^(^)/(i) - hA(s)s. (104) 
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In the limit L — > oo, the effective time evolution operator for Fujikawa's formulation |l9[] is 

BMsr}) ^ expi-^ f F(s„s,)dt\. (105) 



Next, let us examine the adiabatic change of s. Let f{s) be an adiabatic basis for H(s). Now 

H'^is) = f(syH(s)f(s) (106) 

is a diagonal matrix whose non-zero elements are the eigenvalues of H{s). Thanks to the adia- 
batic theorem, we employ the diagonal approximation for Bd({if)): 

BMsr}) - B,Ms,}) = expj--^ j| F^'is, s)dtj , (107) 

where F^(s, s) is defined as 

F^(s, s) = H°(s) - M°(s)s. (108) 
Hence £ad({*r)) is decomposed into geometric and dynamical factors: 

Bad({*,)) = expj/ J A°(i)c!!ijexpJ-i^ //°(s)flf?j . (109) 

Now it is trivial to apply our formulation explained in the main text to obtain the holonomy 
matrix (Eq. dZSTl). 

B. A derivation of Eq. (6Si 



To evaluate Eq. mzi . we need to obtain the image of Ce in the sphere (0, tp). First, we remind 
that is a periodic function of from Eqs. (149 al l and (|59] l. Hence @\lZo must be a multiple of 
In. Second, we remind that there is no spectrum degeneracy along the path Cg if r+ i Z, where 
r± = {A(2 - p)+ ii{2 - q)] I {2k). Namely, the nondegenerate regions are divided into squares by 
the lattice (r+,r_) e 1?. Within each square, 0|g![o takes a constant value. Accordingly, it is 
suffice to evaluate them at representative points of the squares. 

Let us examine the case (r+-r_)/2 is an integer A:, which implies A = 2 cos"'((-l)*cos(7r(r+-i- 
r_)/2),cos0 = (-l)*cos6'sin(;r(r+-Hr_)/2)/sin(A/2), andsin0 = sin6»sin(7r(r++r_)/2)/ sin(A/2). 
Accordingly we have 

e'0 = (_i)[M'-++r_)/2] expliX-l)*^^), (110) 

which implies the image of Cg is also a great meridian loop in the sphere (0, (f). Thus we have 

0|2^o = (-l)*27r. (Ill) 

This result is also valid for the case that ([r+] - [r_])/2 is an integer k. 

Let us examine the case (r+H-r_)/2is an integer/, which implies A - 2cos"'((-l)'cos(;7r(r+- 
r_)/2), cos = (_i)'+K''+-''-)/2] and sin = 0. Namely, the image of Cg is a point in (0, ip) space. 
Accordingly we have 

©Is^o^O- (112) 
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This result is also valid for the case that [r+] - [r_] is an odd integer 
To summarize the argument above, we have, from Eq. (l62t . 



M(Cfl) = exp{/7r(l + [r+]-[r_])). 
which is equivalent with Eq. ( l65T l. 



(113) 



C. Algebraic properties of f „ 

We summarize the algebraic properties of fa. Details are found in, for example, Ref. 
It is straightforward to show that they form a Clifford algebra dTSl l. which impUes the following 
formulas: 



For Ba € ■ 



\\Bf 



(114) 



(115) 



where is a norm of real vector, i.e. ||B|| = Furthermore, for a real unit vector n that 

satisfies ||n|| = 1, 



exp -iA ^ 



- COS A - i 



sin A, 



(116) 



which is 2:7r-periodic in A. Next formula is useful to investigate quantum maps: 

exp (-iAfa) fp exp {-iAfa) 

= (1 - Sap)fp + 5ap{co&{2A)fp - i sin(2/l)). 



(117) 



Finally, we show a representation of fa by complex 4x4 matrices Ta in terms of complex 2x2 
matrices defined in Eq. (l37Ti. i.e.. f„ - f(\T„f^. where fn s ^\et). \Kei), \e2), lA'ei)]: 



To = 



h 

-/2 

T3 







iCTy 

-icTy 




T2 = 










/2 


/2 



, T4 = 




/cr- 


—icr. 








(118) 



D. A derivation of Eq. 



We derive Eq. dTSl l from Eq. (TTSI i. First, let = //(2 - ^)/2 and Bi = A(2 - p)/2. Hence we 



have 



IJ — g-'(j"?+-i/')/2g-!B;,fo/2g-;B,i£„«„f„g-;B^fo/2 



e 



'i{liq+Ap)l2 



(119) 
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where we used Eq. ( II 171 ) above, and, 

k = cos Bjj cos Bi - no sin sin Z?^, 
[() = sin cos + «o cos sin B,i, 
L(#0) s n„ sin 5,1. 

Hence we arrive Eq. ( TTSl ). 



(120) 



E. Diagonalization of "aT^a 

Let «„ (ff = 0, ... ,4) be real, and, normalized, i.e. Yjt=Q - 1- obtain the eigen- 

vectors of 



T({na}) = ^ riata. 



(121) 



Q'=0 



with the help of the quaternionic structure of the Hilbert space induced by the fermion time 
reversal invariance iflTll . We follow the convention of quaternions explained in Ref IitIi . e.g.. 



-1, ij - —ji = k, jk = -kj - i, and ki 



-ik 



J- 



In the Hilbert space "H of spin-|, we employ a right quaternionic action for \tp) e "H: \if/)i = 



i\4i), and, = K(\i//}), which implies \i^}k 
Accordingly, for zj, Wj e C, we have 



[Kim] 



ItA) = Z\\ei) + w\\Kei) -Hzaka) + W2\Ke2) 

= klXzl + + k2)(Z2 + jwi) 



Z\ + jwi 



where 



(122) 



(123) 



is the standard frame of the quaternionic Hilbert space. Hence we obtain a natural correspon- 
dence between four-dimensional complex vector space and two-dimensional quaternionic vector 
space. This induces a representation of t„ by 2 x 2 quaternionic matrix t^, i.e.. 



where 





-1 







" fq ^ff /q 










-j 






-^ 




j 


- -^ 




1 






-i 




1 






i 



k 
-k 



(124) 



(125) 



We introduce spherical variables 0, 77, x (p to parameterize as no - cos0, ni — 
cosx sin 77 sin 0, «2 - sin;^f sin 77 sin 9 n-i - cos (/? cos 77 sin 9, and 714 = sin (p cos 77 sin 9. In terms of 
2x2 quaternionic matrix, we have 



COS (cos 77 - je'^'*'^^^ sin 77)6 sin 9 

e+'^(cos 77 + sin 77) sin 9 — cos 



(126) 
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where we introduce 

h = je'^'^^^\ 

Since = je'^f'^>^^e~'^f*^^j = - 1, we have e'"' = cos 77 + /i sin rj. Hence we have 



t({n„}) = /q 



e-'"'e-'>sin( 



COS0 

e+iVghn^^^Q -cos 9 



(127) 



(128) 



Now it is straightforward to see the eigenvalues of f{[na)) are +1. Let be corresponding 
eigenvectors: 



(-sinf) 



Instead of the two above, we put a phase factor on them in the following: 



(129a) 
(129b) 

(130) 



where we need to take care about the noncommutativity of multiplication in quaternions. In the 
complex Hilbert space, |^±) are expressed as 



6 6 
|^+) = \di)cos- + \d2)sva-. 



I 6\ 
l^-> = l^^i) I - 2/ '^^^ 2' 

where \d\) and \d2) are orthonormal 

1^/2) = lei) I) + \Ke2} (e"'('^-^>/2 ^j^^ | j ^ 



(131a) 
(131b) 

(132a) 
(132b) 



We explain the rest of eigenvectors of f ({no,}) for the complex Hilbert space. They are obtained 
by the time-reversal operation on |^±): 



\K^+) = \Kdi) cos ^ + \Kd2) sin ^, 
\K^-) = \Kdi) (- cos ^) + \Kd2) cos ^, 



where 



\Kdi) = ki) sin I) + |/s:ei) Je+'<*'+^)/2 cos |) , 

Because of A'^ = -1, {|^±), |^r^±)} is a complete orthogonal system for 'H. 
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(133a) 
(133b) 

(134a) 
(134b) 



To compute gauge connections for {|^±), \K^±}}, it is useful to summarize the basis transfor- 
mation between 



(135a) 



and 

/o=[|ei>, \Kei), kz), \Ke2)]. 
It is straightforward to obtain the following 

/ = /oexp(-^g4^jexp(-^g3'^]exp(-^g2?7] exp(-^gi6lj 
where ^q, are 4 x 4 Hermite matrices 








-ih 




-cr, ' 











cr. 





1 = 


ih 





- gi = 


a-. 


- ^3 = 





-0"; 


, 84 = 





cr. 



It is also useful to express ga in terms of t„: gi 

g4 = -!TlT2. 



-IT0T3, g2 = IT1T3, g3 



(135b) 
(136) 

(137) 
-iTiT4, and 
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